Intermittency in forced two-dimensional turbulence 
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We find strong evidence for intermittency in forced two dimensional (2D) turbulence in a fiowing 
soap film experiment. In the forward enstrophy cascade the structure function scaling exponents 
are nearly indistinguishable from 3D studies. Intermittency corrections are present in the inverse 
energy cascade as well, but weaker. Stretched exponential tails of the velocity difference probability 
distribution functions and shock like events at large velocity differences also resemble 3D studies. For 
decaying turbulence, where only the forward enstrophy cascade remains, all signs of intermittency 
disappear. 

PACS numbers; 47.27.Gs, 68.15.+e 



Energy dissipation in three dimensional (3D) turbu- 
lence is punctuated by intermittent bursts W and it has 
been conjectured that intense vortex filaments are re- 
sponsible for these bursts Q . Since vortex filaments can- 
not exist in two dimensions (2D) it is perhaps not sur- 
prising that recent experiments by Paret and Tabeling 
on the inverse energy and forward enstrophy cascades 
of 2D turbulence have found no signs of intermittency 
Simulations by Smith and Yakhot Q for the in- 
verse energy cascade agree with these experiments, but 
simulations by Babiano, DubruUe, and Frick |^ do not. 
The latter simulations go beyond the experiments in that 
they explore not only the isolated inverse energy and for- 
ward enstrophy cascades but also simultaneous cascades. 
In each case Babiano et al. find intermittency only in the 
energy cascade, which is partly at odds with the exper- 
iments. Further work could resolve such discrepancies. 
In this letter we report on 2D turbulence experiments in 
flowing soap films which can probe an isolated enstrophy 
cascade, as done by Paret and Tabeling Q , or simultane- 
ous cascades as simulated by Babiano et al. The presence 
of intermittency in 2D would call for fundamentally new 
ideas about the source of the phenomenon in 2D and 
perhaps even in 3D. 

Kraichnan |Q| proposed that there are two scaling 
regimes in the energy spectrum of isotropic homogeneous 
2D turbulence. Energy conservation, accompanied by 
the assumption that the energy should depend only on 
the wavenumber k and the energy dissipation rate per 
unit mass e leads to, E{k) = Ce^/^/c^^/^. This fol- 
lows the same dimensional arguments first applied to 
3D turbulence by Kolmogorov ||^. In 2D there is a fur- 
ther constraint. The mean square vorticity, or enstrophy 
(il = 1/2 |V X vp), must also be conserved. Through the 
same considerations as above, namely the dependence on 
k and the enstrophy dissipation rate per unit mass rj, 
we find E{k) = C'rf'/^k~^ . This implies the existence 
of a second cascade, the enstrophy cascade. In order for 
the two cascades to be present simultaneously, Kraichnan 
proposed a forward enstrophy cascade and an inverse en- 
ergy cascade 0. 

It is in these two inertial ranges that we look for the 



presence of intermittency. A particularly useful tool is 
the velocity difference, or increment, between two points 
separated by a vector r, <5t;(x,r) = v(x-l-r)— v(x). In the 
present experiment we take both v and r along the di- 
rection of the mean flow. Using this quantity it becomes 
possible to probe the statistics of the flow as a function of 
length scale r. If the turbulent velocity field is self sim- 
ilar with respect to r, then the probability distribution 
functions P((5w(r)) should scale with r, as was the case 
in the experiments of Paret and Tabeling ||^,^ . Intermit- 
tent events break the self similarity of the flow fleld and 
subsequently lead to P{Sv{r)) which do not scale with r. 

This deviation from self-similarity is also manifest 
in the moments of the velocity differences, collectively 
known as the structure functions, 
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We also define a second function, Gp{r) ~ {\5v{r)\P), 
which is often calculated for odd p when a scaling range 
in Sp{r) is narrow or absent. Though the theoretical 
implications of Gp(r), with p odd, are not yet clear we 
include it for completeness. In any case, results derived 
from Godd('") follow exactly the same trend as those de- 
rived from Sevenir) = Gevenir). 

The only exact result involving the structure functions 
is the 2D equivalent of Kolmogorov's four-fifths law, re- 
lating the third order structure function to the energy 
dissipation rate It gives, S'3(r) = |er. Though no 
analytic relation has been derived for structure functions 
of higher orders, self-similarity implies that the structure 
functions should scale like Sp{r) cx r''^ , with C,p = p/3 in 
the energy cascade. Dimensional analysis also suggests 
that S3 (r) oc and (p = p in the enstrophy cascade . 
However, because we will look at the sca ling exponents 
using the extended self-similarity (ESS) [|0| hypothesis 
in the analysis of the structure functions, the results ex- 
pected in the absence of intermittency will be the same 
in the enstrophy cascade as those in the energy cascade. 

The present experiment makes measurements on a 
fiowing soap film suspended between two taut vertical 
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nylon wires using an apparatus similar to that of previ- 
ous experiments Two combs, each 64 cm long, are 
placed in the film so that they form an inverted wedge 
1.2 cm apart at the top and 8 cm apart at the bottom. 
The cylindrical teeth have an average spacing of 1.6 mm 
and a diameter of 0.22 mm. The geometry of the combs 
and the positions where data was collected are shown on 
the right of Fig. [l[ This particular orientation of the 
combs was found to produce the most isotropic and ho- 
mogeneous turbulence, with a nearly Gaussian velocity 
distribution and a turbulence intensity low enough to rely 
on Taylor's frozen turbulence assumption. 
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FIG. 1. Energy spectra at locations A though D as noted 
at right. The energy and enstrophy cascade ranges are with 
respect to curve A, the fully forced case. Theoretically pre- 
dicted energy and enstrophy scaling exponents are provided 
for comparison with the data. 

A laser velocimetry system (LDV) from TSI Inc. is 
used to make the velocity measurements. Data is taken 
for two hours in the center of the film at each of four ver- 
tical locations to obtain data sets of 33.6 million points 
each. Roughly 14 kilometers of film flows past the mea- 
surement location during the taking of each data set, 
with velocities being recorded every 0.4 mm on average 
(with a typical data rate of 5kHz). The mean velocity 
of the film at location (A) in Fig. |l| is w = 214 cm/s 
with a root mean square deviation (Aw^)^/^ = 21 cm/s 
(where Av = v — v). The turbulence intensity is therefore 
It = 9.8% at point A, less than half of values reported in 
other work on 3D turbulence jlj,|l^. At locations pro- 
gressively farther below the combs we see a correspond- 
ing change in the magnitude of the turbulence intensity, 
decreasing to 8.0% at B, 6.8% at C, and 6.0% at D. 

From the data sets, the P{6v{r)) were calculated at 
various values of r, and from these, the functions Gp{r) 
were calculated by numerical integration akin to Eq. ^. 
It should be noted that there is an inherent limitation on 
the maximum order of the structure function that can be 
calculated in any physical experiment. As p is increased 



the majority of the weight of the integrand comes increas- 
ingly from the tails of the PDFs, where the number of 
data points becomes quite small. We found that the ma- 
jority of the weight of the numerical integration of Gp{r) 
was not contributed by the wings of the PDFs for p < 14. 
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FIG. 2. Normalized velocity difference probability distribu- 
tions for a range of separations through the energy and enstro- 
phy cascades from data taken at location A. The thick gray 
line is a Gaussian reference curve. The upper inset shows one 
particularly steep velocity fluctuation event representative of 
the wings of the distributions. The lower inset is an average 
of all such events with \5v/ {5v'^}^^^\ > 5.35 at a separation of 
0.5 cm (there were a total of 311). 

Furthermore, as suggested by Anselmet et al. the 
convergence of the structure functions with an increase 
in the size of the data set provides a check on the quality 
and quantity of the data. If each small subset of the data 
is statistically identical to all other subsets, we expect to 
arrive at the same value for the Gp{r) regardless of the 
size of the subset used in the calculation or which subset 
in particular is chosen. We find that for all separations 
r, the variation of Ge{r) is below 5% for data sets of 
5 • 10^ points or greater. For Gi2(r), 25 • 10^ points are 
needed to meet the same requirement. Given the above 
considerations we are confident of the measured structure 
functions up to order 12. 

The energy spectrum in Fig. |^ shows simultaneous ev- 
idence of the two scaling regimes, with E{k) cx k^^/^ in 
the inverse energy cascade (associated with driven turbu- 
lence) and E{k) oc k~^ in the forward enstrophy cascade 
(associated with decaying turbulence). The location of 
the bend between the two scaling ranges gives an esti- 
mate of the effective injection length scale at each loca- 
tion. We estimate it to be A;/27r = 1/A = 0.63 cm~^ or 
A = 1.6 cm at location A. Subsequent figures will refer 
to the energy and enstrophy cascades which are inferred 
from the spectrum at A. More details on the turbulent 
spectra and the nature of the injection scale were previ- 
ously pubhshed by Rutgers [ pl| . 
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FIG. 3. The exponent a of the best fit stretched exponen- 
tial function P{Sv) oc exp {~\Sv/a\"} (for data from location 
A). The solid line is the fit for 5v/{5v^)'^^^ > 2 (i.e. right 
wing); the dashed line corresponds to Sv/(6v^y^^ < —2 (left 
wing) . 

With the energy and enstrophy cascade inertial ranges 
defined, we inspect the data for signs of intermittency. 
Recent experiments on the separate energy and enstro- 
phy cascades of a quasi 2D experimental system by Paret 
and Tabeling found Gaussian velocity difference 

PDFs independent of the scale and concluded that there 
was a complete absence of intermittency. From the analy- 
sis of our forced turbulent data (measured at point A) , we 
conclude quite the opposite. Figure ^ shows P{dv{r)) for 
r covering both the forward enstrophy and inverse energy 
cascades. To accentuate relative shape changes the width 
of each distribution was rescaled by its second moment. 
The distributions with r in the energy cascade are close to 
Gaussian and differ substantially from distributions with 
r in the enstrophy cascade, which are decidedly more ex- 
ponential in nature. This variation of the shape can be 
quantified by fittin g th e distributions to a stretched ex- 
ponential function P{Sv) oc exp{-|§^|"}. The 
dependence of the stretching exponent a on r suggests a 
flow-field which is not self-similar and is an indication of 
intermittent turbulent flow. The positive and negative 
tails of the PDFs beyond 2a are then fit independently 
1^,^. The trend toward a slower than Gaussian falloff 
at small length scales is shown in Fig. ||. The minimum 
a for negative Sv is greater than for positive 5v reflect- 
ing the inherent asymmetry. Note that the sign of the 
asymmetry is in keeping with the exact theoretical result, 
S,{r) = ler. 

The events in the raw velocity data which populate 
the wings of the P{6v{r)) can be singled out. A typical 
event is shown in the upper inset of Fig. |2[ The lower 
inset shows an average of all such events for r = 0.5 
cm and \5v{r)/ {Sv{r)^y^'^\ > 5.35 (an arbitrary cutoff). 
Since the third moment of P{Sv{r)) is positive, the av- 
erage has a corresponding characteristic sharp negative 
slope. Note that the insets in Fig. || have been plotted as 



raw time traces from the measurement probe, as is cus- 
tomarily done in previous work. The application of the 
frozen turbulence assumption will lead to a shock with 
positive slope which gives a positive S3 (r) for our defini- 
tion of 6v{r). The shock-like nature of these events bears 
a strong qualitative resemblance to intermittent events 
measured in 3D turbulence, but typically with a velocity 
increment of opposite sign [^jl6), since S3{r) is negative 
for 3D turbulence. 




2 4 6 8 10 12 
P 



FIG. 4. A. ESS scaling exponents in the inverse energy (A) 
and forward enstrophy (Q) cascades at location (A) (see Fig. 

The solid line is the K41 prediction, and the dashed line 
the log-Poisson. Error bars show the 95% confidence interval. 
B. Measurements of Cp,ESS in enstrophy range at points 
A,B,C, and D (see Fig. jlf 

These intermittent events have a strong effect on the 
structure functions Sp{r), which theoretically scale with 
r. Practically, it is extremely difficult to discern this scal- 
ing directly since one typically requires several decades 
of clearly discernible scaling in E{k). Most experiments, 
ours included, do not have that luxury. The standard 
remedy follows a technique introduced by Benzi et al. 
]lOt where Sp{r) is plotted against G3(r), leading to much 
cleaner power laws throughout the inertial range. This 
procedure is referred to as extended self-similarity and is 
used in most studies of intermittency in 2D and 3D tur- 
bulence. We will refer to the scaling exponents derived 
from this technique as Cp.ESS- Figure ^A shows the Cp,ess 
as a function of the order p, for data collected at point A. 
The solid line is the prediction in the absence of intermit- 
tency, and the dashed curve the log-Poisson prediction 
[]l7| (which is a good fit to the canonical 3D intermit- 
tency measurements of Anselmet et al. Our scaling 
exponents from the enstrophy cascade of forced 2D tur- 
bulence (O) are nearly indistinguishable from those for 
3D measurements. Exponents from the energy cascade 
deviate as well from the self-similar prediction, but to a 
lesser degree (A). 

At first sight our measurements appear to contra- 
dict those of Paret and Tabeling mentioned earlier . 
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There is, however, an important difference between our 
measurements and theirs. Paret and TabeUng's electro- 
convection ceh is hmited by its range of accessible length 
scales and is capable of producing only the inverse energy 
cascade or the forward enstrophy cascade independently. 
Computer simulations are often similar in scope when 
they employ hyperviscous terms to suppress the enstro- 
phy range when the energy range is being studied 
The results shown in Fig. ^ are for a system with both 
cascades present simultaneously. 

To verify the importance of two cascades versus one, we 
performed the analysis leading to Fig. ||A again for points 
progressively further downstream from the combs where 
the turbulence has decayed (points B, C, and D in Fig. 
|l|). The results are summarized in Fig. IJb. At the largest 
distance below the comb (D) the spectrum indicates only 
an enstrophy cascade ||ll|,|l^,|l^ . Accordingly the indica- 
tions of intcrmittency disappear with Cp.ess tending to- 
ward p/3. Our results at this location thus agree with 
the electroconvection experiments Q on the isolated en- 
strophy cascade. We were not able to produce the case 
of an isolated energy cascade in the flowing soap films 
and thus, cannot make a comparison to the correspond- 
ing electroconvection experiments |^ which also found 
no intcrmittency Soap film and electroconvection exper- 
iments therefore agree for similar situations, but both dis- 
agree with parts of simulations by Babiano et al, which 
claim intcrmittency for inverse energy cascades but never 
for forward enstrophy cascades, regardless of whether the 
cascades are isolated or simultaneous. 

There is always a question of the viability of flowing 
soap films as effective 2D incompressible fluids. It is pos- 
sible that our results are influenced by quantities specific 
to flowing soap films, such as air drag, fllm thickness 
variations, or compressibility effects. Though these ef- 
fects are important and we will make them the subject 
of further study, the clear identification of the energy and 
enstrophy cascades in flowing films jll|,|2^ justifies this 
experimental medium. To this can be added our current 
agreement with electroconvection experiments that 
there is no significant intcrmittency in 2D turbulent flows 
with an isolated enstrophy cascade. 

In conclusion, find evidence of intcrmittency in forced 
2D turbulence when both the inverse energy and forward 
enstrophy cascades are present simultaneously. Signs of 
intcrmittency are strongest in the enstrophy cascade, and 
present to a lesser degree in the energy cascade. Evidence 
of intcrmittency is threefold. First we find a considerable 
dependence of the shape of the velocity difference PDFs 
on the scale r. As r is increased through the enstrophy 
and energy cascade ranges the tails of P{6v{r)) can be fit 
by a stretched exponential with exponent a varying from 
1.23 toward Gaussian statistics. Secondly we can identify 
shock like changes in the velocity which contribute to the 
tails of P{Sv{r)). Thirdly, there is a deviation of the ESS 
structure function scaling exponents from the K41 pre- 



diction in both the energy and enstrophy cascades. The 
experimental values lie well-below the K41 prediction of 
Cp — P/3 in the energy cascade, and in the enstrophy 
cascade lie directly on top of the log-Poisson prediction 
jlTf . The same system shows none of these signs of intcr- 
mittency when the turbulence is left to decay, in which 
case there is only an isolated enstrophy cascade. This 
observation agrees with other experiments using an elec- 
tronconvection technique Q . The simultaneous presence 
of both inverse energy and forward enstrophy cascades 
appears necessary to observe intcrmittency in 2D turbu- 
lence. 
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